Dimensionality effects in the LDOS of ferromagnetic hosts probed via STM: 
spin-polarized quantum beats and spin filtering 
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We theoretically investigate the local density of states (LDOS) probed by a STM tip of ferro- 
magnetic metals hosting a single adatom and a subsurface impurity. We model the system via 
the two-impurity Anderson Hamiltonian. By using the equation of motion with the relevant Green 
functions, we derive analytical expressions for the LDOS of two host types: a surface and a quantum 
wire. The LDOS reveals Friedel-like oscillations and Fano interference as a function of the STM 
tip position. These oscillations strongly depend on the host dimension. Interestingly, we find that 
the spin-dependent Fermi wave numbers of the hosts give rise to spin-polarized quantum beats in 
the LDOS. While the LDOS for the metallic surface shows a damped beating pattern, it exhibits 
an opposite behavior in the quantum wire. Due to this absence of damping, the wire operates as a 
spatially resolved spin filter with a high efficiency. 

PACS numbers: 07.79.Fc, 74.55. +v, 85.75.-d, 72.25.-b 



I. INTRODUCTION 

The local density of states (LDOS) of electronic sys- 
tems with impurities can exhibit Fano line shapes, due 
to the quantum interference between different electron 
paths. Such interference arises from the itinerant elec- 
trons that travel through the host and tunnel into the 
impurity siteSfii 2 - For a single magnetic adatom in the 
Kondo regime 3 - probed by a scanning tunneling micro- 
scope (STM) tip, interesting features manifest when one 
has a spin-polarized electron bath present. Here we men- 
tion the splitting of the Kondo peak in the differen- 
tial conductance due to the itinerant magnetism of the 
host4Such hallmark has already been found experimen- 
tally in an Fe island with a Co adatom^ Additionally, the 
STM system can also operate as a Fano- Kondo spin-filter 
due to a spin-polarized tip and a nonmagnetic hosti&iln 
the absence of a ferromagnetic host, the Fano- Kondo 
profile becomes doubly degenerated— Away from the 
Kondo regime, a spin diode emerges^ 

In the condensed matter literature on scanning mi- 
croscopy, there is a profusion of work discussing spin- 
dependent phenomena employing ferromagnetic leads 
coupled to quantum dots or adatoms in the Kondo 
regime i 4 ' 6 ' 7 ' 31 " ^ Here we mention those with metallic 
samples and buried impurities, in which the anisotropy 
of the Fermi surface plays an important role in electron 
tunneling According to the experiment of Priiser et 
al2^, such anisotropy allows atoms of Fe and Co beneath 
the Cu(lOO) surface to scatter electrons in preferential di- 
rections of the material due to an effect called "electron 
focusing". In this scenario, the STM becomes a new tool 
for the detection of the Fermi surface signatures in the 
real lattice of a metal. In contrast, much less attention 



has been devoted to spin-polarized systems away from 
the Kondo regime and with two impurities. 

Thus we present in this work a theoretical description 
of the systems sketched in Fig. [TJ We show that in- 
teresting phenomena such as the spin-polarized quantum 
beats in the LDOS and the spin-filtering effect arise. To 
this end, we consider two distinct geometries consistent 
with recent experiments: a metallic surface and a quan- 
tum wire. The 2D case emulates the Fe island in Ref.[f|. 
The quantum wire on the other hand, mimics the "elec- 
tron focusing" effect investigated in Ref. [25]. Interest- 
ingly, we note that the pioneering quantum wire treat- 
ment for "electron focusing" in a side-coupled geometry 
can be found in Ref . [29l| . In this treatment the nonin- 
teracting single impurity Anderson model** was solved 
in one dimension by considering the impurity above the 
wire. We should also point out that the full ab-initio cal- 
culation that yields to "electron focusing" in Ref. [25| can 
be qualitatively recovered by the simple quantum wire 
model adopted in Ref. [29| . 

Here we extend this one-dimensional treatment of the 
Anderson Hamiltonian by including a spin-dependent 
DOS for the wire, a second lateral impurity right be- 
neath it and Coulomb interaction in both impurities. We 
perform our study in the framework of the two-impurity 
Anderson model by employing the equation of motion 
approach to calculate the LDOS of the system. The 
Hubbard I approximation^ 2 - is used by assuming for the 
sake of simplicity infinite Coulomb energies at the impu- 
rities. We show that the LDOS can be written in terms 
of the Fano factor, the Friedel-like function for charge 
oscillations and the spin-dependent Fermi wave numbers 
of the host. Such quantities lead to spin-polarized quan- 
tum beats in the LDOS. We also show that this effect 
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is strongly correlated to the host dimensionality. Thus 
the quantum beats in the LDOS of the metallic sur- 
face present a long-range damped behavior in contrast 
to the undamped one found in the quantum wire system. 
Such distinct features originate from the specific forms 
assumed by the Fano factor and Friedel function, which 
depend on the dimensionality of the host. Therefore the 
metallic surface and the quantum wire become spatially 
resolved spin filters, where the latter displays a higher 
efficiency due to the undamped LDOS. 
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Figure 1. (Color online) Side-coupled geometry with two im- 
purities in presence of a STM tip. Ttip is the tip-host coupling, 
a) left panel: 2D evanescent waves appear in the LDOS of a 
metallic surface. The system is treated as a two-dimensional 
electron gas showed in the right panel, b) left panel: the 
confinement of ID waves in specific directions (perpendicu- 
lar wave fronts) is due to the "electron focusing" effect (see 
Ref.[2^|). Each direction is modeled by a quantum wire as 
illustrated in the right panel. 

This paper is organized as follows. In Sec. [TT1 we show 
the theoretical model of the ferromagnetic hosts with the 
impurities in the side-coupled geometry as sketched in 
Fig. [T]and derive the LDOS formula for both systems, the 
metallic surface and the quantum wire. The decoupling 
scheme Hubbard 1^ for the Green functions is presented 
in Sec. Mil In Sec. IIV1 we discuss the results for the 
quantum beats in the LDOS and the spin-filtering. The 
conclusions appear in Sec. [V] 



II. THEORETICAL MODEL 

A. Hamiltonian 

In order to probe the LDOS of the ferromagnetic hosts, 
we represent a STM tip weakly connected to hosts hy- 
bridized to a pair of side-coupled impurities as outlined 
in Fig. Q] The systems we investigate are described ac- 
cording to the two- impurity Anderson model given by the 
Hamiltonian^ 
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The spin-polarized electron gas forming the hosts is de- 
scribed by the operator ct ( C £ CT ) f° r the creation (an- 
nihilation) of an electron in a quantum state labeled by 
the wave vector k, spin a and energy s? . For the im- 
purities, dj a (dj a ) creates (annihilates) an electron with 
spin a in the state Ejda, with j = 1,2. The third term of 
Eq. {1} accounts for the on-site Coulomb interaction Uj 
at the jth impurity placed at position Rj . In our calcu- 
lations, we assume U± = U2 — > 00 (single occupancy of 
the impurities). Finally, the last two terms mix the host 
continuum of states and the levels Sjd a ■ This hybridiza- 
tion occurs at the impurity sites Rj via the host-impurity 
couplings V-j: a and the plane waves <fi^ a (Rj^ = e** -i ^ . 

A/" CT is the number of conduction states for a given spin 
a. The ferromagnetic hosts are considered spin-polarized 
electron baths, characterized by the polarization 



in which 
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(3) 



is the density of states of the hosts in a Stoner-like 
framewor k 53 i 54 expressed in terms of the spin-dependent 
half-width D a and the density po for the case P = 0. 



B. LDOS for the spin-polarized systems 

To obtain the host LDOS we introduce the retarded 
Green function in the time coordinate, 



x (n 



where 



(4) 



(5) 



is the fermionic operator describing the quantum state of 
the host site placed below the STM tip, h is the Planck 
constant divided by 27r, 9 (t) the step function at the in- 
stant t, j3 = 1/ksT with ks as the Boltzmann constant 
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and T the system temperature, Zfm an d \ n ) are the par- 
tition function and a many-body eigenstate of the system 
Hamiltonian [Eq. (JTJ)], respectively, and [• • ■ ,■••]+ is the 
anticommutator for Eq. ([5]) evaluated at distinct times. 
From Eq. ([4J the spin-dependent LDOS at a site R of the 
host [see Fig. [I] can be obtained as 



(6) 



where Q a (e + , R) is the time Fourier transform of G a (t, R). 
Here, e + = e + ir\ and r/ — > + . In what follows we first 
develop a general formalism for impurities localized at 
arbitrary positions Rj and Ri ; later on we take the limit 
Rj = Ri = 0, in order to treat the side-coupled geometry 
of this work. 

To obtain an analytical expression for the LDOS we 
apply the equation-of- motion approach to Eq. ((3]). Thus 
we substitute Eq. (0 in Eq. ((4]) and begin the procedure 
with 

*> 3) = TT E <^ (*) ^ (*) (7) 
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Performing ^ on Eq 
d 
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where we have used 
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In the energy coordinate, we solve Eq. © for 
GZ~r* an d obtain 
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(11) 



Notice that we need to find the mixed Green function, 
Gd c -( £+ )- To this end, we define the advanced Green 
function 



{n\[4 a {0),c&(tj\ \n), (12) 



which results in 



(n\ [d]^ (0) , c 9 - CT (i)] |n) - -e^^ (t) 
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where we used once again Eq. (|10[) , interchanging fc O <f. 
Thus the Fourier transform of Eq. (|13|) becomes 



^ jC? (e") = ea^e, (e") + ^ E (^) 

x (O > (1 4 ) 

with e~ = e — zt7- Applying the property 

0W e+ ) = {^<* (*")}' (15) 
to Eq. (fl4|) . we show that 
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and 
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Now we substitute Eq. (fTT|) into Eq. (fTTj) and the 
latter into Eq. JTJ) in the energy coordinate to obtain 
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It is worth mentioning that the imaginary part of the first 
term of Eq. (fT5|) gives the background DOS of the host 
[Eq. (J3j)] and the others describe impurity contributions, 
with 
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being the Fano parameter due to the single coupling 
between the host and a given impurity. This factor en- 
codes the quantum interference originated by electrons 
traveling through the ferromagnetic conduction band 
that tunnel to the impurity state and return to the band, 
and those that do not perform such trajectory. Addition- 
ally, we recognize 



A >° = - V E V fkA fa) Ha fa 6 {* ~ 
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as an expression that we call the Friedel function, be- 
cause it leads to Friedel-like oscillations in the LDOS with 
OLjcr as a spin-dependent phase. At the end, the Green 
function Q a (e + ,R) indeed depends on G% d ( £ ) anc ^ the 
mixed Green function G^. d (e). Finally, from Eqs. (J5]) 
and (fT5)h the LDOS of the ferromagnetic systems can be 
recast as the expression 
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where 

&Qjlo = -gjaQUr^^djd, ( £ )} + i\ A ja\qia 

x cos (a ja + 7^) + \ A ja\ \Ai a \ COS (Ctja - a Zcr ) 
- q ja \A Xa \ cos (a !(T + |j Imj^d, 0)} 



+ fea \ A ia\ sin (a; CT + -J + qia \A ja \ 

x sin (ctja + tj) + |A,-| sin (a,- CT - a icr ) 

x Re{^. d! (e)}. 
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The set of Eqs. ([2T|) and (|22[) is the main analyti- 
cal finding of this work. It describes the spin-dependent 
LDOS in ferromagnetic hosts with two impurities local- 
ized at distinct sites i?, . In the absence of the last term 



of Eq. (|2"Tj) , it reduces to the case of two decoupled sys- 
tems with one impurity each. The terms in Eq. (|22p. 
indeed, hybridize such single- impurity problems via the 
mixed Green functions G d a ( £ )- Thus the LDOS for- 
mula encodes the single Fano factor qj a . the Friedel-like 
function Aj a and the new interfering term Agjia- 

We close this section by recalling that the phase a.j a 
is nonzero for the quantum wire system, as we shall see 
later. The quantities qj a and Aj a in the quantum wire 
device are indeed functions that exhibit undamped oscil- 
lations as the tip moves away from the impurities. Con- 
versely, damped oscillations are predicted in the metallic 
surface setup. In this case OLj a = and Aj a becomes a 
real function. Thus the quantity \Aj a \ should be read 
just as Aj a in Eqs. (f2"Tj) and (f2"2"|) . Moreover, a ferromag- 
netic environment is characterized by two spin-dependent 
Fermi wave numbers, namely, fcp-j- and kFi, which at low 
polarization P introduces a slightly difference between 
them. As a result, this feature leads to a full LDOS 
Pldos Pldos with spin-polarized quantum beats, that 
can be damped or undamped, depending upon the sys- 
tem dimensionality. We shall look more closely at these 
features later. 

In STM experiments, in particular within the linear 
response regime and neglecting tip-adatom coupling, the 
differential conductance G = + is the observable 
measured by the tip, whose spin component is given by^ 



G a = — ttF 
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where e is the electron charge (e > 0), Tu v is the tip- host 
coupling, / is the Fermi-Dirac distribution and <p is the 
applied bias. For <fi < the host is the source of electrons 
and the tip is the drain. For <j> > 0, we have the opposite. 
It is useful to define the dimensionless LDOS 



LDOS 



t 4. 
Pldos j~ Pldos 



Pfm\ + Pfmi 
and the transport polarization, 



(24) 



Pt = 



G 1 ' - G l 
G^ + G^' 



(25) 



in order to investigate the spin-polarized quantum beats 
as we shall see in Sec IIV1 Recall that, in the absence 
of the impurities, the transport polarization of Eq. (|25j) 
becomes Pt = P, as established by Eq. ([2]). In Sec. 
Wl we shall verify that Eq. (f2"5|) oscillates around P, 
exhibiting two distinct behaviors as a result of the system 
dimensionality - damped spin-polarized quantum beats 
in the metallic surface setup and an undamped pattern 
in the quantum wire device. 
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C. Fano and Friedel-like functions for the metallic 
surface system 

In this section, we calculate the expressions for the 
Fano parameter [Eq. (fL9|) ] and the Friedel-like func- 
tion [Eq. ([20]) ] in the metallic surface case, where no 
manifestation of "electron focusing" occurs. This calcu- 
lation was previously performed in the single-impurity 
problem^ and now we present an extension applied to 
the double impurity system. We begin by solving Eq. 
(|20| . To this end, we assume ^(R) = e ikRcosek " for the 
electronic 2D wave function and use 
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the angular representation for the zeroth-order Bessel 
function. Thus, in the wide-band limit \e\ -C D a with 
the flat-band DOS 
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k—kf 



expressed in terms of the spin-dependent Fermi wave 
number kp a and an element of area S in the host, we 
find 



(32) 



which allow us to close the calculation. Equations (|3"Tj) 
and (]32[) imply the relation 
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As the Friedel function is already known from Eq. ([28JI , 
the quantity —^Gja (e, R) provides a relationship for the 
Fano parameter. To this end, we can write 



1 = 1 PFMcr 

7T/90 2 po 
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with Gjia {s, R) obeying the following integral represen- 
tation: 
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In the equation above we have used, for the sake of 
simplicity, the linear dispersion relation 



£/cct — D a kp a {k kFa) 



(36) 



A 3a = P -^L\ 
Po 



(kpoRj 



= A 2D 



R — Ri 



(28) 



as the rel- 



for the Friedel-like function with R = 

ative coordinate with respect to the jth impurity. Notice 
that, according to Eq. (|20[) . phase ctj a is zero and Eq. 
(|28|) is a real quantity. In the case of the Fano parameter, 
we start defining the advanced Green function 



k 



17] 



by assuming the Lorentzian shape 



R 



(29) 



the Hankel functions if^ (£) = Jo(£) + *^o(0 ancl 

H o 2 \0 = MO - iYo(£)- We remark that Eq. (JSgJ in 
combination with Eq. (0) for the ferromagnetic host al- 
lows us to find 



l + P 



(37) 



In particular, for a small polarization P, Eq. (]3T|) re- 
sults in slightly different Fermi wave numbers and, con- 
sequently, in spin-polarized quantum beats in the full 
LDOS as we shall see. Looking at Eq. |35|) . we calculate 
the integral Gjia by choosing a counterclockwise contour 
over a semicircle in the upper-half of the complex plane, 
which includes the simple pole e^ a = +iA. Applying the 
residue theorem, we obtain 
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for an energy-dependent coupling^ in order to obtain an 
analytical solution for qj a . Notice that, in the limit A 3> 
\ska\, Eq. |30|) recovers the case V-^ a = V. Thus we can 
write the identities 



with &a = kp a (l + i-^-)- For the evaluation of Gj2& we 
used a clockwise contour over a semicircle in the lower- 
half plane, including the poles e^a = £ — ir) and e^a = 
—iA, which yields 
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and 
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with k F 



k Fa (l 



erty (£) = (C 



, Eq. 



becomes 



Taking into account the prop- 
for the second term in Eq. 



with m as the effective electron mass and £ as a given 
length in the wire. Additionally, in the wide-band limit 
|e| <C D a , we find the following complex Friedel-like func- 
tion: 
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Explicit calculation of the terms in the brackets of 
Eq. (gDD leads to 



A, 
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kpa-R. 



characterized by a spin-dependent phase a 
This phase results in nondamped, oscillatory behavior as 
a function of R. which also appears in qj^ ■ Thus we take 
into account Eqs. (|4*4l and (|4"5"j) . rewriting Eq. (fTT)|) as 
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where we have assumed \e\ <C D a , A <C D a and A ^ |e|. 
In order to ensure the limit V^ a = V in Eq. (|3"0)) , we 
make the substitution of Eqs. fll|, (gDJ), (gT|) and (gU) 
in Eq. (|3"3")l , showing that 



with 
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in the limit |e| <C D a and with p as the principal value. 
Finally, we obtain the Fano factor 



C = o 
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ID 



for any value of k Fa R- 

In summary, the zero value of the Fano parameter 
given by Eq. (|43|) and the zeroth-order Bessel function 

Jo (k Fa -R\ found in Eq. (|28|) lead to long-range damped 

spin-polarized quantum beats in the full LDOS. This fea- 
ture will be discussed in Sec. IIVI 



D. Fano and Friedel-like functions for the quantum 
wire system 

Here we determine the Fano parameter in Eq. ([TO]) and 
the Friedel-like function in Eq. (|20[) for the quantum wire 

case. Following A. Weismann^ we use 4>k<j(R) = e%kR as 
the electron wave function, in which the direction intro- 
duced by R defines the STM tip-impurity distance where 
"electron focusing" manifests. We also use the dispersion 
relation 
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h 2 k 2 
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which also presents spin-dependent Fermi wave numbers 
k F -f and k F ± as in Eq. ([4"6"| . Here they are still connected 
via Eq. (|37p . thus leading to undamped spin-polarized 
quantum beats in the full LDOS. 



III. CALCULATION OF THE IMPURITY 
GREEN FUNCTION 



In the present section we calculate G dd[ (e) (j, I = 1,2) 
that appear in Eqs. flU]) and (25} for the LDOS. To 
handle the interacting term of the Hamiltonian, we adopt 
the Hubbard I approximation^ which provides reliable 
results at temperatures above the Kondo temperature^ 
We begin by repeating the equation-of-motion approach 
for these Green functions, which results in 



(e - e jda + iTj ja )Q^. d . = 1 + UiQd jo v. ijS ,d ja 

+£( E «*- a v) ><su ( 5 °) 
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where Sj d(T = £jda + ^fja f° r ' 7^ 3- I n the equation 
above, Gd la n dllil d ja is a higher-order Green function ob- 
tained from the time Fourier transform of 



x ( n | (t)n«j, ff (t),dt ff (0) |n), (52) 

L J + 

with ti^ct = dj^dig being the number operator of the Ith 
impurity with spin a (opposite to a). Here 



E 



Ri 



(53) 



and apply the equation-of-motion approach to 
C ,t , . . In order to cancel the second term 

with the last one in Eq. (|55[) , we combine the approxi- 
mations in Eqs. (|57j) and (|58[) simultaneously with the 
property 



E W = E zmz v kM-o W , (59) 



which is only fulfilled on a metallic surface system, while 
for the quantum wire, it holds in the side-coupled geom- 
etry Ri = Rj = 0. Bearing this in mind, we can rewrite 
Eq. (|55f as follows 



represents the real part of the noninteracting self-energy 
and 



(e - eua -Ui + irj) Gd la n dlS ,d ja = Sij {n dlS ) 



^ = = -^E^%A~, (4) c 0) + f E r^sA* (#) I x ^^-(eo) 



5 ( £ - £ fc J 



(54) 



describes the corresponding imaginary part, which plays 
the role of a generalized Anderson parameter In 
order to close the system of Green functions in Eqs. (|50|) 
and (|5"Tj) . we first take the time derivative of Eq. (|5^|) 
and then perform the time Fourier transform. With that 
we obtain 



(e + - side, - Ui) Qd^n dlS ,d ja = $lj {n dl s) 

+ (e (^)) x Ka**»*. 

(55) 

which also depends on new Green functions of the same 
order of Gdi„n d s ,d ja an( ^ on * ne average occupation num- 
ber (n dl g), that is, calculated as 

(»*»>= /" +00 &(--Imfe d! )l/( £ ). (56) 



Within the Hubbard I approximation we truncate the 
reen functions G i 

Cldlgdlvtdja 

ing to the decoupling scheme 



Green functions G t and G,t , , accord- 

^cldi s di„,dj„ 3 d\ s c%-di a ,dj a 



c\_d ls di^,d j 



g dle^A.-\ C k* dl *) g U 



cy i9 )Gl d] (57) 



(58) 



Once again, employing the equation-of-motion approach 
for G jt j j , we find 



^ktr) ^cr d^-die ,djcr 



•E r ^-^r x C (%) aw*,- (6i) 

Here we continue with the Hubbard I scheme, proceed- 
ing as in Eqs. (|5"T|) and ([55)) by making the following 
approximations: 



dls C q<j) Gc z d j<7 , 



$c\ s d ls c % „,d ia ~ \4 S C q<y)Gc %a d D ^ 
Qd- ja n dlS ,d ja - {ndis) Gd-.dj 

and replacing Eq. (|55|) in Eq. (|6"Tj) to show that 



(62) 
(63) 
(64) 



Gdin d s ,dj 



A; it 7 



e+ -s r . 



x (n dl a)Gl d .. 



(65) 
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To close the original setup of Green functions in Eqs. 
([50]) and ([ST]), we substitute Eq. (JB5J) in Eq. (J5UJ and 
obtain 

+ (n^) ^ (sg, (e) - iT- jUT ) x C^ d ., (66) 



where 



rL = ^J2\V rka \ 2 5(e-s^), (67) 



which allows us to determine all the necessary Green 
functions for the LDOS. Thus, to solve the system com- 
posed by Eqs. ([50]) . ([5"T]) and ([66]) . we now assume the 
side-coupled geometry Ri = Rj = (Fig. [JJ). A ge- 
ometry with Ri 7^ Rj will be published elsewhere. 
By choosing the side-coupled configuration and assum- 
ing constant symmetric couplings V^ a = V^ a = V, we 
verify from Eqs. (J5J) and ([57]) that r, JCT = T lja . = T° Ua = 
T a = r p ™ g depends on the standard Anderson param- 
eter F = irV 2 po- Additionally, in the wide-band limit, 
Eq. ([53]) ensures T,f ja = S§ CT =0. For the sake of sim- 
plicity we consider the infinite Coulomb correlation limit 
(U\ = U2 — > 00). Thus the direct Green function for the 
impurity j = 1 reduces to the form 



Gt ldl (e) 



1 - {ndis) 



£ ~ t\2do + iAi2 



(68) 



IV. NUMERICAL RESULTS 

A. Numerical Parameters 

Here we present the results obtained via the formu- 
lation developed in the previous section. The energy 
scale adopted is the Anderson parameter V. We em- 
ploy the following set of model parameters: T = 0.2 eV, 
£ida = £id = -lOr and e 2da = e 2d = -4.5TJ&2± Such 
values correspond to a Kondo temperature Tk ~ 50-ff 
found in the system Co/Cu(lll) with Coulomb interac- 
tion U = 2.9 eV^^- Thus the Hubbard I approxi- 
mation is employed with T = r/lOfce = 231. IK just to 
avoid Kondo physics. Finally, in order to generate spin- 
polarized quantum beats in the LDOS, we substitute Eq. 
© with P = 0.1 in Eq. (57J). 



— 1 1 1 1 

metallic surface 




where 



sir 



£\1do — £\da + Si2cr (E) 



(69) 



1.1 



represents a renormalized energy level dressed by the real 
part of the nondiagonal self-energy 



and 



Si2<r (e) = - (1 - (n dlS )) (1 - (n d2 a)) 

(g ~ £2da) p 2 

(£~£2d*) 2 + ri ° 



Ai 2ff = T a - (1 - (n dl a)) (1 - (n d2 a)) T c 



(70) 



(£ - £2da) + Tl 



(71) 



is an effective hybridization function. The mixed Green 
function G^ 2( j 1 (e) becomes 

ft*M--M ( .'^i) }ft*M- (72) 

Notice that the other Green functions G% ^ and Q% d 
can be derived by swapping 1 O 2 in Eqs. $38|) and ([72}. 
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Figure 2. (Color online) In both panels we use ksT = O.ir. 
(a) LDOS [Eq. ([2i]l] of a metallic surface with P = 0.1 as a 
function of e/Y for different values of IzfR in the short-range 
limit (see panel (b)). The Fano profile presents two antireso- 
nances placed at e = ei<j = — 10r and e = £2d = — 4.5r, which 
display an evanescent behavior for increasing distances, (b) 
Keeping the energy at e = — 10T, Friedel oscillations appear 
in the LDOS. 
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Figure 3. (Color online) In both panels we use ksT = O.ir. 
(a) LDOS [Eq. (TJ2])] of a metallic surface with P = 0.1 as a 
function of e/Y for different values of UfR in the long-range 
limit (see panel (b)). The Fano profile presents two antireso- 
nances placed at e = Eid = — 10r and e = Eid = — 4.5r, which 
display an oscillatory behavior for increasing distances, (b) 
Damped spin-polarized quantum beats emerge in the LDOS 
as function of kpR with e = Eid = — 10F. 
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Figure 4. (Color online) (a) LDOS [Eq. {24])] at k B T = O.ir 
of a quantum wire with P = 0.1 as a function of e/T for dif- 
ferent values of positions kpR- Fano-profiles appear around 
e — eid = — 10F and e = Eid = — 4.5r. This pair of antireso- 
nances (black line) can be tuned to resonances (green line) by 
moving laterally the STM tip. (b) In opposition to the metal- 
lic surface device, undamped spin-polarized quantum beats in 
the LDOS at e — Eid = — 10F manifest. 



B. Metallic Surface 

Defining kpi = kp, we begin the analysis in the metal- 
lic surface apparatus by dividing our study into regions 
we call short, intermediate and long-range limits. The 
short-range limit presented in Fig. Uta) reveals that the 
LDOS given by Eq. (|24p as function of energy exhibits 
two Fano antiresonances. Each one corresponds to the 
discrete levels of the adatom (j = 1) and the subsurface 
impurity (j = 2). The main feature in this situation is 
that the Fano profile conserves its line shape when the 
dimensionless parameter kpR is changed. Additionally, 
this profile is suppressed for increasing distances, tending 
to the DOS background of the host. 

In Fig. [DJb) we look at how the LDOS evolves with 
kpR exactly at the Fano antiresonance e = = — 10r. 
At the host site (fcpi? = 0), the LDOS presents a deple- 



tion. Such a dip in the LDOS is a result of charge screen- 
ing around the impurities by conduction electrons, which 
suppresses the LDOS of the host. Beyond the adatom 
position, the LDOS is indeed dictated by Friedel oscilla- 
tions, which also lead to a strong decay in the long-range 
limit [see the Fig. (Ub)]. The evanescent feature of the 
LDOS is a result of the interplay between the Friedel- 
like expression A 2 ^ and the Fano parameter qfP . These 
quantities are governed by Eqs. (|28p and (|43 )) . where 
the former evolves spatially according to the zeroth-order 
Bessel function Jo- Such damping in the LDOS has been 
already observed experimentally in a system composed 
by an Fe host and a Co adatom^ 

In Fig. [3la) we plot the Fano line shape in the long- 
range regime (kpR > 15). The same dips at e = — 10T 
and e = — 4.5r are observed as in the short-range limit 
[Fig. Ufa)]. However, a contrasting feature is found be- 
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tween these two limits. While in the short-range case the 
dips become suppressed as kpR increases, in the long- 
range limit the dip oscillates with kpR. This is a result 
of the oscillatory profile observed in the LDOS for in- 
creasing kpR [see Fig. [U^b)]. The oscillations of the dip 
can be more clearly visualized in Fig. E^b), where we 
show the LDOS at e = £<ji = — 10r. A peculiar beat- 
ing is observed in the LDOS due to the slightly different 
Fermi wave numbers [see Eq. (|37l) ]. 

C. Quantum Wire 

Figure HJa) shows the LDOS plotted against energy 
for different k F R values in the short-range limit. For the 
STM tip at k F R = 0, the LDOS shows the two-dip struc- 
ture already observed in Fig. [U^a). In contrast, as kpR 
increases, the antiresonances change to resonances, pass- 
ing through intermediate profiles (asymmetric Fano line 
shapes). We emphasize that this behavior in the LDOS 
was recently observed in the experiment performed by 
Pruser et at with atoms of Fe and Co beneath the 
Cu(100) surfaced 

We observe in Fig. 0|b) the evolution of the LDOS 
with kpR. Non-evanescent oscillations occur, modulated 
by an amplitude beating. This undamped behavior is 
encoded by Eqs. PS|) and for Friedel-like oscil- 

lations (Ajp ) and Fano interference (qj^), respectively. 
These quantities are simple trigonometric functions with- 
out damping. This feature is due to the absence of an 
extra dimension for the scattering of the electronic wave. 
On the other hand, in 2D this propagation is spread in a 
plane leading to a spatial decay in the LDOS. Thus the 
amplitude of the undamped beats is much larger than in 
the metallic surface device. This means that in such case, 
the LDOS signal can be more easily resolved experimen- 
tally. 



D. Transport Polarization 




20 40 60 80 100 120 140 



~T~\ — ' — I — 1 — I — ' — I — 1 — I — ' — I — 1 

(b) quantum wire 




.0.6 I ■ I ■ I ■ I ■ I ■ I ■ I ■ 



20 40 60 80 100 120 140 

k p R 

Figure 5. (Color online) Transport polarization of the ferro- 
magnetic hosts [Eq. (|25[) ] at fcsT = O.ir and applied bias 
</> = eid — — 10T. (a) Damped spin-polarized quantum beats 
appear in the polarization of the metallic surface device, (b) 
Undamped type occur in the polarization of the quantum 
wire. In both situations, we have a spatially resolved spin- 
filter, with a polarization that oscillates around P = 0.1. In 
the quantum wire case, this oscillation is more pronounced. 



Another quantity we investigate is the spin- 
polarization given by Eq. (|2"o]l . As the differential 
conductance of Eq. (|2"5| is proportional to the LDOS 
[Eq. (pM|) ]. the ferromagnetic hosts filter electrons that 
tunnels into (or out of) the STM tip. This filtering 
is dominated by the majority spin component. Thus 
devices without impurities behave as spin filters with a 
spatially uniform polarization that coincides with the 
value given by Eq. Here we adopt P = 0.1. Due 

to the impurities in the side-coupled geometry and the 
host dimensionality, this polarization is perturbed in 
two different forms. In both the metallic surface and the 
quantum wire as we can see in Figs. 03a) and [S](b) with 
applied bias = E\d = — 10r, the polarization oscillates 
around Pt = P = 0.1. Unlike the metallic surface 
system, where small deviations with damping occur, the 
"electron focusing" effect in the quantum wire leads to 



undamped and pronounced oscillations. The polariza- 
tion in the latter case does not exceed Pt ~ +0.62 or 
fall below Pp ~ —0.5. Therefore the polarized current 
through the junction formed by the STM tip and the 
surface alternates from spins up (+0.62) to down (—0.5) 
depending on the tip position. Additionally, along 
this probing direction, the polarization not only can 
invert the orientation of the majority spin component, 
but also becomes zero at some sites, where locally the 
unbalance of spins is totally suppressed. As a result, we 
have a tunneling current without polarization in specific 
positions on the sample surface. On the other hand, as 
we can see in Fig. Eta), the amplitude of the beats in the 
metallic surface polarization is extremely suppressed and 
does not change its signal (Pt > 0). Thus the quantum 
wire operates as a spatially resolved spin filter, with a 
higher efficiency. 
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V. CONCLUSIONS 



In order to investigate a ferromagnetic system with 
two impurities, we have calculated, the LDOS and the 
spin polarization of hosts in two different dimensional- 
ities. Impurities in the side-coupled geometry, as out- 
lined in Fig. [TJ were taken into account. We analyzed 
both a metallic surface and a quantum wire described 
by the two-impurity Anderson model in the picture of 
a spin-polarized electron gas with impurities away from 
the Kondo regime. We presented a model in which an un- 
perturbed ID electron host in the presence of localized 
states produces undamped behavior in the LDOS Fano 
profile [see Figs. HJa) and@|b)], similar to that observed 
experimentally^ In contrast, our 2D model revealed a 
damped oscillatory behavior [Figs. [UJa), [UJb), and 
[31(b)] . We demonstrated that these opposed features orig- 
inate from the interplay between the Friedel-like function 
and the Fano parameter, which assume different func- 
tional forms according to the host dimensionality. Keep- 
ing the energy fixed and tuning the STM tip position, we 



verified the emergence of spin-polarized quantum beats 
in the LDOS given by Eq. (|2~4")) as well as in the transport 
polarization of Eq. ([25]) . Such an effect is due to interfer- 
ence between the slightly different Fermi wave numbers 
kpf and kp± [Eq. (|5Tj) [ in the LDOS, achievable in hosts 
with low spin polarizations. Therefore the quantum wire 
setup behaves as a spatially resolved spin-filter with a 
high efficiency, as we can see in Fig. [SJb). Away from 
the adatom, this device can magnify or invert locally the 
original spin orientation of the host, also displaying sites 
where this polarization is completely quenched. As a pos- 
sible experimental implementation of this apparatus, we 
suggest the systems investigated by Pruser^ Such se- 
tups present the same one-dimensional character as our 
effective quantum wire model. 
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